The Effective Stability of Adaptive Timestepping ODE Solvers

نویسنده

  • Harbir Lamba
چکیده

Abstract. We consider the behaviour of certain adaptive timestepping methods, based upon embedded explicit Runge-Kutta pairs, when applied to dissipative ODEs. It has been observed numerically that the standard local error controls can impart desirable stability properties, but this has only been rigorously verified for very special, low-order, Runge-Kutta pairs. The rooted-tree expansion of a certain quadratic form, central to the stability theory of RungeKutta methods, is derived. This, together with key assumptions on the sequence of accepted timesteps and the local error estimate, provides a general explanation for the observed stability of such algorithms on dissipative problems. Under these assumptions, which are expected to hold for ‘typical’ numerical trajectories, two different results are proved. Firstly, for a large class of embedded Runge-Kutta pairs of order (1, 2), controlled on an error-per-unit-step basis, all such numerical trajectories will eventually enter a particular bounded set. This occurs for sufficiently small tolerances independent of the initial conditions. Secondly, for pairs of arbitrary orders (p−1, p), operating under either error-per-step or error-per-unit-step control, similar results are obtained when an additional structural assumption (that should be valid for many cases of interest) is imposed on the dissipative vector field. Numerical results support both the analysis and the assumptions made.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Dynamical Systems and Adaptive Timestepping Ode Solvers

Initial value problems for ODEs are often solved numerically using adap-tive timestepping algorithms. We formulate a large class of such algorithms as discrete dynamical systems which are discontinuous and of higher dimension than the underlying ODE. By assuming suuciently strong nite-time convergence results on some neighbourhood of an attractor of the ODE we prove existence and upper semicont...

متن کامل

Stepsize control for the Milstein scheme using first-exit-times

We introduce a variable timestepping algorithm for the Milstein scheme applied to SDEs with a scalar stochastic forcing. Multiple local error estimates are used, corresponding to different terms appearing in the Taylor series of the local truncation error. The timesteps are then chosen so as to bound the standard deviation of the contribution of each of these terms over a unit time interval, an...

متن کامل

Stabilisation of Discrete Adjoint Solvers through Improved Primal Timestepping

One major challenge in applying the adjoint method to large industrial cases is the lack of robustness of the flow and adjoint solvers. Discrete adjoints transpose the Jacobian of the system of discretised equations exactly, and hence inherit the linear stability from the primal flow solver. Typically the linear stability properties of the iterative scheme of the flow solver is very well unders...

متن کامل

Stability and convergence of a finite volume method for two systems of reaction-diffusion equations in electro-cardiology

The monodomain equations model the propagation of the action potential in the human heart : a very sharp pulse propagating at a high speed, which computation require fine unstructured 3D meshes. It is a non linear parabolic PDE of reaction diffusion type, coupled to one or several ODE, with multiple time-scales. Numerical difficulties, such as unstructured meshes and stability are addressed her...

متن کامل

An Instance of Failure for the MATLAB Explicit ODE45 Solver

We consider the adaptive strategies applicable to a simple model describing the phase lock of two coupled oscillators. This model has been used to show an instance of failure of the ODE45 RungeKutta-Felberg solver implemented within the MATLAB ODE suite, see [J. D. Skufca. Analysis still matters: a surprising instance of failure of Runge-KuttaFelberg ODE solvers. SIAM Review, 46:729-737, 2004]....

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Numerical Analysis

دوره 43  شماره 

صفحات  -

تاریخ انتشار 2005